J. Jpn. Bot. 

70 : 322 - 327 ( 1995 ) 
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Levels of genetic diversity and population structure in Hosta venusta and H. tsushimensis , herbaceous 
perennials endemic to Cheju Island and Tsushima Island, respectively, were investigated using allozyme 
markers. Although the distributions of the species are restricted in the islands, they maintain high levels of 
genetic variation (mean He = 0.296 and 0.215 for the populations of H. venusta and H. tsushimensis, 
respectively) and the degrees of population differentiation were lower than the species of Hosta with similar life 
history traits and widespread distribution range such as H. capitata and H. clausa. The findings suggest that H. 
venusta and H. tsushimensis may have large effective population sizes. 


Recent reviews of plant allozyme literature have 
shown causal relationships between life history and 
ecological traits and the amount and patterns of ge¬ 
netic variation (Hamrick and Godt 1989, Hamrick, 
Godt and Sherman-Broyles 1992). For example, spe¬ 
cies with endemic distributions maintain lower level 
of genetic diversity than species with widespread 
distributions. However, several studies reported that 
geographic range is not always a good predictor of the 
genetic structure of species (Karron et al. 1988, 
Olmstead 1990, Soltis and Soltis 1991, Godt and 
Hamrick 1993, Chung and Chung 1994). This may be 
largely explained by two possible explanations. First, 
species, even congeners, often differ in many other 
aspects of their biology. Second, relatively fewer 
allozyme studies of endemic or rare species, particu¬ 
larly island endemic species, are available in the plant 
allozyme literature. 

In this study I investigate levels and distribution of 
genetic variation in Hosta venusta F. Maekawa and H. 


tsushimensis N. Fujita, herbaceous perennials re¬ 
stricted to Cheju Island, Korea and Tsushima Island, 
Japan. Hosta venusta occurs commonly on grassland 
and understory of reforestation areas, while H. 
tsushimensis grows on diverse habitats such as forest 
margin, valleys and rocky areas (Fujita 1976; Chung 
pers. obs). Like other hostas, these two species are 
insect-pollinated. Flowers are visited by bees and 
seeds with wing are dispersed by wind. Chromosome 
number is 2n = 60 (Chung et al. 1991a). 

Materials and Methods 

A total of 124 samples of H. venusta and H. 
tsushimensis rootstocks were collected from four lo¬ 
calities in Cheju and Tsushima Islands. Table 1 lists 
the four populations, their alphabetic codes, sample 
sizes, and vouchers used in this study. Collections of 
rootstocks from each site were made at least 2 m apart 
to reduce the likelihood of collecting clones. The 
rootstocks were transported to the Botany Plant Growth 
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Table 1. Population code, location of the four populations of//. 
venusta and H. tsushimensis 


Code 

Locality 

N a 

H. venusta 

A 

Cheju Island, Cheju city, Korea 

32 

B 

Cheju Island, Samgumburi, 
Pukcheju-gun, Korea 

31 

H. tsushimensis 

C 

Tsushima Island, Izuhara, Pref., Japan 

27 

D 

Tsushima Island, Mt. Ariake, Japan 

34 


“Sample size. 


Facility at the University of Georgia where plants 
were grown under uniform conditions for electro¬ 
phoretic analyses. Voucher specimens of all collec¬ 
tions were deposited at GA, GNUC, KYO, MO, and 
SNU. 

Leaf tissue for each sample was crushed in liquid 
nitrogen with a mortar and pestle, and proteins were 
extracted using a phosphate polyvinylpyrrolidone 
buffer (Mitton et al. 1979). Protein extracts were 
absorbed onto Whatman 3 mm chromatography pa¬ 
per wicks, placed in 96-well microtest plates and 
stored at-70°C. Electrophoresis was performed using 
11% starch gels. Except as noted, gel and electrode 
buffers and enzyme staining procedures from Soltis et 
al. (1983) were used to assay ten enzyme systems: 
isocitrate dehydrogenase (IDH) on system 2; shikimate 
dehydrogenase (SKDH) on system 4; phosphogluco- 
mutase (PGM), diaphorase (DIA), and glutamate 
oxaloacetate transaminase (GOT) on system 7; leu¬ 
cine aminopeptide (LAP), triosephosphate isomerase 
(TPI), /Lgalactosidase (/3 GAL), and fluorescent este¬ 
rase (FE) on a modification (Haufler 1985) of system 
8. The staining procedures for DIA and /3 GAL 
followed the method described by Cheliak and Pitel 
(1984). Putative loci were designated sequentially, 
with the most anodally migrating isozyme designated 
1, the next 2, and so on. Likewise, alleles were 
designated sequentially with the most anodally mi¬ 
grating alleles designated a. 


Levels of allozyme variation were estimated at the 
species level and for individual populations with a 
computer program developed by M. D. Loveless and 
A. Schnabel. A locus was considered polymorphic if 
two or more alleles were observed. Genetic diversity 
was measured by four parameters: percent polymor¬ 
phic loci (P), average number of alleles per locus (A), 
mean effect number of alleles per locus (Ae), and gene 
diversity or expected heterozygosity (He). The statis¬ 
tics of the four parameters were well described in 
Hamrick and Godt (1989) and Hamrick, Godt, and 
Sherman-Broyles (1992). 

Observed heterozygosity,' expected heterozygos¬ 
ity, and Wright’s (1922) fixation indices ( F ) were 
calculated for each locus and for each population to 
reveal significant deviations in genotype frequencies 
from expectations under Hardy-Weinberg equilib¬ 
rium. Positive values of F indicate an excess of 
homozygotes and negative values indicate an excess 
of heterozygotes within the populations. Significance 
of the F coefficients was tested by a % 2 -statistic 
following Li and Horvitz (1953). 

Allozyme variation was partitioned within and 
among populations for each locus using Nei’s (1973, 
1977) genetic diversity formulae. The total gene di¬ 
versity for polymorphic loci (H T ) was partitioned into 
within population (H s ), and among population (D sx ) 
components. The proportion of genetic diversity found 
among populations was calculated as G sx = D sx /H x . 
G sx values were tested for significance by % 2 = 2AG SX 
(a-1); df = (a- 1) (n- 1), where A is the total number of 
individuals analyzed, a is the number of alleles at the 
locus and n is the number of populations (Workman 
andNiswander 1970). Nei’s (1972) genetic identities 
and distances were calculated for pairwise combina¬ 
tion of populations and species. 

Results 

Sixteen and 11 of 24 loci were polymorphic in H. 
venusta and H. tsushimensis , respectively. Three alleles 
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were observed for five loci of the two species, with the 
total average numbers of alleles for the two species 
being 1.88 and 1.42, the average effective number of 
alleles 1.57 and 1.40, and overall genetic diversity 
0.315 and 0.208. Populations of H. venusta and H. 
tsushimensis maintained high levels of genetic varia¬ 
tion: the mean percentage of polymorphic loci of 66 
and 45% and the mean genetic diversity of 0.296 and 
0.215 (Table 2). 


Fifteen (eight positive; seven negative) of 32 fixa¬ 
tion indices in H. venusta and 11 (seven positive; four 
negative) in H. tsushimensis were significantly differ¬ 
ent from zero (Table 3). On the whole, allele frequen¬ 
cies were close to Hardy-Weinberg expectations with 
mean F IS of -0.034 and 0.027 for H. venusta and H. 
tsushimensis. 

Mean total genetic diversity at the polymorphic 
loci (H t ) for H. venusta (0.473) was very comparable 


Table 2. Summary of allozyme variation for 24 loci within populations of H. venusta and H. 
tsushimensis'. proportion of polymorphic loci (P), mean number of alleles per locus 
(A), mean effective number of alleles per locus (Ae), mean observed heterozygosity 
(Ho), and mean expected heterozygosity (He). 


Code 

P 

A 

Ae 

Ho (SD) 

He (SD) 

H. venusta 

A 

66.67 

1.83 

1.57 

0.272 (0.065) 

0.286 (0.048) 

B 

66.67 

1.88 

1.58 

0.344 (0.064) 

0.305 (0.046) 

Mean 

66.67 

1.86 

1.57 

0.308 (0.046) 

0.296 (0.033) 

H. tsushimensis 

C 

45.83 

1.42 

1.44 

0.193 (0.056) 

0.215(0.049) 

D 

45.83 

1.42 

1.43 

0.186 (0.051) 

0.216(0.048) 

Mean 

45.83 

1.42 

1.44 

0.189 (0.038) 

0.215 (0.034) 


Table 3. 

Fixation indices (F) for 16 polymorphic loci of H. venusta and 11 polymorphic loci of H. 


tsushimensis 





H. venusta 

H. tsushimensis 

Locus 





A 

B 

C 

D 


Idh 

-0.4318** 

-0.3864* 

-0.5588*** 

-0.5581*** 

Skdh-2 

-0.1250 ns 

0.5770** 

0.3690 ns 

—0.1917 ns 

Pgm-1 

-0.0862 ns 

-0.2449 ns 

- 

- 

Pgm-2 

-0.2353 ns 

-0.2000 ns 

-0.3590 ns 

-0.5227*** 

Pgm-3 

-0.3404 ns 

-0.2449 ns 

0.4692** 

0.6036*** 

Dia-1 

0.3848* 

0.1712 ns 

- 

- 

Dia-2 

-0.5000** 

-0.6053*** 

-0.0192 ns 

-0.3958* 

Dia-3 

-0.1982 ns 

-0.6138 ns 

0.7017*** 

0.4591** 

Got-1 

- 

- 

0.1452 ns 

0.5550*** 

Got-2 

0.5714*** 

0.5841*** 

—0.145 l ns 

0.0337 ns 

Got-3 

0.0762 ns 

-0.6620*** 

-0.3701 ns 

-0.0828 ns 

Pgd-1 

0.4337** 

-0.2931** 

0.0063 ns 

0.1434 ns 

Lap-3 

0.5333*** 

-0.1200 ns 

- 

- 

Tpi-2 

-0.002 P s 

0.6155*** 

0.4960*** 

0.6240*** 

Tpi-3 

0.1703 ns 

-0.2380 ns 

- 

- 

Fe-2 

0.6921*** 

0.6013*** 

- 

- 

Fe-3 

-0.1564 ns 

-0.5250*** 

- 

- 
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Table 4. Nei’s (1973, 1977) statistics of genetic diversity for 16 polymorphic loci in H. venusta and 11 polymorphic loci in H. 
tsushimensis. A % 2 test for allele frequency heterogeneity between populations is given: * = P < 0.05; ** = P < 0.01; 
*** = />< 0 . 001 . 


Locus 



H. venusta 



H. tsushimensis 


H t 

H s 

D st 

G st 

x 2 

H t 

H s 

D st 

G S t 

x 2 

Idh 

0.4213 

0.4210 

0.0002 

0.0006 

0.074(2) ns 

0.4656 

0.4656 

0.0000 

0.0000 

0.001(2) ns 

Skdh-2 

0.4274 

0.3571 

0.0703 

0.1645 

20.722(1)*** 

0.4020 

0.3970 

0.0051 

0.0126 

1.537(l) ns 

Pgm-1 

0.2561 

0.2494 

0.0067 

0.0262 

3.305(1 ) ns 






Pgm-2 

0.4058 

0.3987 

0.0071 

0.0174 

4.863(2) ns 

0.4350 

0.4322 

0.0028 

0.0064 

0.782(2) ns 

Pgm-3 

0.3628 

0.3612 

0.0016 

0.0043 

0.543(3) ns 

0.4192 

0.4159 

0.0033 

0.0079 

1.592(3) ns 

Dia-1 

0.4717 

0.4424 

0.0293 

0.0620 

7.816(1)** 






Dia-2 

0.4636 

0.4627 

0.0009 

0.0020 

0.255(1 ) ns 

0.2956 

0.2630 

0.0260 

0.1103 

13.458(1)*** 

Dia-3 

0.4898 

0.4897 

0.0001 

0.0002 

0.024(2) ns 

0.4988 

0.3708 

0.1280 

0.2566 

31.299(2)*** 

Got-1 






0.5575 

0.5525 

0.0050 

0.0089 

1.740(2) ns 

Got-2 

0.4919 

0.4076 

0.0844 

0.1715 

21.605(2)*** 

0.4020 

0.3973 

0.0047 

0.0118 

1.439(2) ns 

Got-3 

0.5200 

0.5074 

0.0126 

0.0242 

4.321 (2) ns 

0.3651 

0.3425 

0.0226 

0.0620 

9.866(2)** 

Pgd-1 

0.6269 

0.6089 

0.0180 

0.0287 

8.447(2)* 

0.5169 

0.5068 

0.0101 

0.0195 

5.271 (2) ns 

Lap-3 

0.4874 

0.3431 

0.1443 

0.2961 

37.306(2)*** 






Tpi-2 

0.4819 

0.4776 

0.0043 

0.0088 

1.112(2) ns 

0.6202 

0.6201 

0.0001 

0.0001 

0.066(2) ns 

Tpi-3 

0.5668 

0.5600 

0.0068 

0.0119 

2.294(2) ns 






Fe-2 

0.4717 

0.4496 

0.0221 

0.0468 

5.899(2) ns 






Fe-3 

0.6159 

0.5565 

0.0594 

0.0964 

32.992(2)*** 






Mean 

0.4726 

0.4433 

0.0292 

0.0601 


0.4525 

0.4331 

0.0195 

0.0451 



with that for H. tsushimensis (0.453) (Table 4). The 
proportions of allozyme diversity among populations 
(G S t) of H. venusta and H. tsushimensis were 0.060 
and 0.045, indicating low levels of genetic differen¬ 
tiation among populations. Only six of 16 polymor¬ 
phic loci for H. venusta and three of 11 polymorphic 
loci for H. tsushimensis revealed significant differ¬ 
ences in allele frequencies among populations (Table 
4). Average genetic identities among the H. venusta 
and//, tsushimensis populations were 0.934 and 0.974. 
Mean genetic identities among the two species was 
0.764. 

Discussion 

The review of allozyme diversity in plant 
populations by Hamrick and Godt (1989) provides 
information for interpreting the allozyme data pre¬ 
sented here. Hosta venusta and H. tsushimensis ex¬ 
hibit considerably more genetic variation than that for 
plant species with similar life history taits (e.g., long- 


lived herbaceous species with a predominantly ani¬ 
mal-outcrossing mode of reproduction and wind- 
dispersed seeds, and endemic geographic range). Mean 
percent polymorphic loci within populations of spe¬ 
cies with these life history traits is 35%, mean number 
of alleles per locus is 1.53, and 1.15 mean effective 
number of alleles per locus. In addition, the mean 
values of genetic diversity (He) within populations of 
H. venusta and H. tsushimensis (0.296 and 0.215) 
were substantially higher than those for species with 
similar life history traits (0.113). The similar results 
were observed in H. yingeri S. B. Jones. Although this 
species is only restricted to remote islands off the 
southwestern coast of Korea such as Taehuksan, 
Sohuksan, and Hong Islands (Chung and Kim 1991), 
it maintains higher levels of genetic variation within 
populations (He = 0.250; Chung and Chung 1994). 

In general, geographic range has been shown to be 
significantly associated with the levels of variation 
maintained within populations (Hamrick, Linhart and 
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Mitton 1979, Hamrick and Godt 1989). Widely dis¬ 
tributed plant species tend to maintain more variation 
than more narrowly distributed species. Contrary to 
the expectation, H. venusta and//, tsushimensis, spe¬ 
cies with a restricted geographic range, have higher 
levels of genetic variation than mean estimates of 
genetic diversity within populations for more wide¬ 
spread their congeners such as H. capitata (Koidz.) 
Nakai (He = 0.153; Chung 1994a) and H. clausa 
Nakai (He = 0.082; Chung 1994b). Species with 
small, isolated populations might be expected to have 
less genetic diversity than species with large and 
contiguous populations (Chung and Kang 1994). Most 
populations of H. capitata and H. clausa are small 
and isolated (Chung et al. 1991a), which seems to be 
one of the major factors contributing to lower levels 
of genetic variation within the species. 

The relatively high levels of genetic variation 
found in H. venusta and H. tsushimensis is consistent 
with several aspects of their biology. Hosta venusta 
and H. tsushimensis appear to be predominantly 
outcrossing because flowers had spatial separation 
(herkogamy) of mature anthers and stigma impeding 
self-pollination. Hosta venusta and H. tsushimensis 
are relatively common in Cheju andTsushimalslands 
(Fujita 1976; Chung pers. obs.). It has been observed 
that seeds are dispersed by wind and pollinator visita¬ 
tion areas are found to be large, which consistent with 
the results that, on the whole, allele frequencies in 
populations of both species were close to Hardy- 
Weinberg expectations. Species that evolve from 
hybridization or through polyploidy should have en¬ 
hanced levels of genetic variation. The chromosome 
number is relatively high (2n = 60; Chung et al. 
1991a), H. venusta and H. tsushimensis were sug¬ 
gested as a polyploid origin (Chung et al. 1991a). The 
combination of these factors may contribute to main¬ 
tain high levels of genetic variation within the two 
species. It is of interest to note that populations of H. 
minor (Baker) Nakai, a Korean endemic species oc¬ 


curred in the eastern and southern Korean Peninsula, 
maintain similar levels of genetic variation with those 
of H. venusta and H. tsushimensis (He = 0.230; Chung 
1994c). This may support that the three species are of 
similar phylogenetic histories (Chung 1991, Chung et 
al. 1991 b), or ecological and life history traits (Chung 
1994c). 

Genetic differentiation among populations is prin¬ 
cipally a function of gene flow among populations via 
pollen and seeds dispersal (Ellstrand and Marshall 
1985). Of the total variation in H. venusta (G sx = 
0.060) and H. tsushimensis (G ST = 0.045), 6% and 
slightly less than 5% are due to difference among 
populations. This is supported by relatively homog¬ 
enous allele frequencies between populations for the 
two species. Only six and three out of 16 and 11 
polymorphic loci for H. venusta and H. tsushimensis, 
respectively, were found to be significantly different 
in allele frequencies between populations. The levels 
of divergence found in H. venusta and H. tsushimensis 
are substantially lower than the means for species 
with endemic geographic range (G sx = 0.248), 
outcrossing-animal breeding system (G sx = 0.197), 
seed dispersal by wind (G sx = 0.143), and long-lived 
herbaceous perennials (G sx = 0.213) (Hamrick and 
Godt 1989). However, too few populations were 
surveyed for each species to draw any conclusions 
about estimating gene flow by indirect method. 

In summary, it is highly probable that the higher 
levels of variability and lower degrees of differentia¬ 
tion found within and between populations of H. 
venusta and H. tsushimensis than in species with 
similar life history traits would suggest that higher 
effective population sizes are maintained in the spe¬ 
cies resulting from extensive gene flow via seed 
dispersal by wind, large population sizes, and high 
fecundity (Chung and Chung 1994). 

I thank J. L. Hamrick, S. B. Jones, W. Anderson, 
R. Wyatt, S. Sherman-Broyles, N. Fujita, and S. S. 



December 1995 


Journal of Japanese Botany Yol. 70 No. 6 


327 


Kang for assistance. Special thanks go to K. Kaneko 
for his advice on field trips in Tsushima Island. I am 
also grateful to the staffs of the Department of Na¬ 
tional Monuments of Culture Property Preservation 
Bureau, Korea for permissions to collect samples 
from Mt. Halla, Cheju Island, Korea. This research 
was supported in part by an U.S. National Science 
Foundation Dissertation Improvement Grant (BSR- 
8914430) and a Korea Research Foundation, non 
directed research fund (1992). 

References 

Cheliak W. M. and Pitel J, A. 1984. Techniques for starch gel 
electrophoresis of enzymes from forest species. Informa¬ 
tion Report PI-X-42, Petawawa National Forestry Institute, 
Agriculture Canada, Canadian Forestry Service, pp. 1-49. 
Chalk River, Ontario. 

Chung M. G. 1991. Notes on phylogenetic relationships among 
species of Korean Hosta (Liliaceae) based on chloroplast 
DNA restriction data. Plant Species Biol. 6: 71-74. 

- 1994a. Genetic structure in Korea populations of 

Hosta capitata (Liliaceae). J. Plant Biol. 7: 277-284. 

- 1994b. Low levels of genetic diversity within 

populations of Hosta clausa (Liliaceae). Plant Species Biol. 
9:177-182. 

-1994c. Genetic variation and population structure in 

Korean endemic species: III. Hosta minor (Liliaceae). J. 
Plant Res. 107: 377-383. 

-and Chung H. G. 1994. Allozyme diversity and popu¬ 
lation genetic structure in Korean endemic species: II. Hosta 
yingeri (Liliaceae). J. Plant Biol. 37: 141-149. 

——-, Hamrick J. L., Jones S. B. and Derda G. S. 1991a. 

Isozyme variation within and among populations of Hosta 
(Liliaceae) in Korea. Syst. Bot. 16: 667-684. 

--, Jones S. B., Hamrick J. L. and Chung H. G. 1991b. 

Morphometric and isozyme analysis of the genus Hosta 
(Liliaceae) in Korea. Plant Species Biol. 6: 55-69. 

——-and Kang S. S. 1994. Genetic variation and population 

structure in populations of Eurya japonica (Theaceae) in 
Korea. Amer. J. Bot. 81: 1077-1082. 

-and Kim J. W. 1991. The genus Hosta Tratt. (Liliaceae) 

in Korea. Sida 14: 411^120. 

EllstrandN. C. and Marshall D. L. 1985. Interpopulational gene 
flow by pollen in wild radish, Raphanus sativus. Am. Nat. 

r > '> Hosta venusta i 7 2 7 4 4 0 is H. tsushi- 
mensis 

^ lt ©&mmm 


126: 606-616. 

Fujita N. 1976. The genus Hosta (Liliaceae) in Japan. Acta 
Phytotax. Geobot. 27: 66-96 (in Japanese). 

Godt M. J. W. and Hamrick J. L. 1993. Genetic diversity and 
populations structure in Tradescantia hirsuticaulis 
(Commelinaceae). Amer. J. Bot. 80: 959-966. 

Hamrick J. L. and Godt M. J. W. 1989. Allozyme diversity in 
plant species. In: Brown A. H. D., Clegg M. T., Kahler A. L. 
and Weir B. S. (eds.), Plant Population Genetics, Breeding 
and Genetic Resources, pp. 43-63, Sinauer, Sunderland. 

-,-and Sherman-Broyles S. L. 1992. Factors 

influencing levels of genetic diversity in woody plant spe¬ 
cies. New Forests 6: 95-124. 

--, Linhart Y. B. and Mitton J. B. 1979. Relationships 

between life history characteristics and electrophoretically 
detectable genetic variation in plants. Ann. Rev. Ecol. Syst. 
10: 173-200. 

Haufler C. H. 1985. Enzyme variability and modes of evolution 
in Bommeria (Pteridaceae). Syst. Bot. 10: 92-104. 

Karron J. D., Linhart Y. B., Chaulk C. A. and Robertson C. A. 
1988. The genetic structure of populations of geographi¬ 
cally restricted and widespread species of Astragalus 
(Fabaceae). Amer. J. Bot. 75: 1114-1119. 

Li C. C. and Horvitz D. G. 1953. Some methods of estimating 
the inbreeding coefficient. Amer. J. Human Genet. 5: 107— 
117. 

Mitton J. B., Linhart Y. B., Sturgeon K. B. and Hamrick J. L. 
1979. Allozyme polymorphisms detected in mature needle 
tissue of ponderosa pine. J. Hered. 70: 86-89. 

Nei M. 1972. Genetic distance between populations. Am. Nat. 
106: 283-292. 

- 1973. Analysis of gene diversity in subdivided 

populations. Proc. Natl. Acad. Sci. USA 70: 3321-3323. 

- 1977. /^-statistics and analysis of gene diversity in 

subdivided populations. Ann. Human Genet. 41: 225-233. 

Olmstead R. G. 1990. Biological and historical factors influenc¬ 
ing genetic diversity in the Scutellaria angustifolia 
(Liliaceae). Evolution 44: 54-70. 

Soltis D. E., Haufler C. H., Darrow D. C. and Gastony G. J. 
1983. Starch gel electrophoresis of ferns: A compilation of 
grinding buffers, and staining schedules. Amer. Fern J. 73: 
9-27. 

Soltis P. S. and Soltis D. E. 1991. Genetic variation in endemic 
and widespread plant species: Examples from Saxifragaceae 
and Polystichum (Dryopteridaceae). Aliso 13: 214-223. 

Workman P. L. and Niswander J. D. 1970. Population studies on 
southwestern Indian tribes. II. Local genetic differentiation 
in the Papago. Amer. J. Human Genet. 22: 24—49. 

Wright S. 1922. Coefficients of inbreeding and relationship. 
Am. Nat. 56: 330-338. 


4- 44 4 '> H. Capitata 4 H. clausa 4 4" 4 ’0 4 4 

tl h . 



